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NECESSARY THICKNESS OF RANDOMLY ORIENTED ARMOR 
FOR METEOROID PROTECTION 
byC. D. Miller 
Lewis Research Center 

SUMMARY 

Weighting factors developed earlier were used to obtain a histogram of velocity and 
a revised influx rate for meteoroids relative to a gravity-free Earth, fully corrected for 
photographic biasing effects of meteoroid mass, real velocity relative to Earth's atmos- 
phere, gravitational focusing, and direction in space of approach of meteor toward Earth. 
Similar results were derived for conditions near the real Earth. 

Bimodal log-normal equations were derived to match the histograms of velocity for 
both gravity- free Earth and real Earth conditions. A double integration yielded an equa- 
tion for general use to determine the necessary magnitude of a damage criterion, such as 
armor thickness, for a desired probability of a successful mission from the meteoroid 
damage standpoint. The double integration involved the bimodal log-normal velocity dis- 
tributions and a well confirmed equation for influx rate of meteoroids of various masses. 
As an example of the method of use of the equation for necessary magnitude of a damage 
criterion, use was also made of an expression for penetration by hypervelocity impact 
provided in a personal communication by Richard H. Fish and James L. Summers of 
Ames Research Center. However, if preferred by the reader, other criteria could 
readily be used in the same manner. 


INTRODUCTION 

Meteoroids, small particles within the solar system and in orbit about the sun, will 
continually impinge on radiators and other components of space vehicles. For space 
missions of long duration, meteoroids having masses within the so-called photographic 
range may represent a serious hazard. It is therefore desirable to develop accurate 
analytical expressions for calculation of probability of meteoroid impact damage. 

Until recently, necessary refinements in analysis of meteoroid influx data have been 



lacking. Some recent efforts directed toward the necessary refinements will be sum- 
marized here, and the manner in which practical use may be made of the results will be 
illustrated. 

Photographic meteoroids, on impact with the upper atmosphere, produce luminous 
trails, or meteors, of sufficient brilliance that they may be photographed. In the pho- 
tography of such meteors, two important biases exist. Particles having greater velocity 
relative to Earth's atmosphere produce trails more easily photographed than those mov- 
ing more slowly. Also, heavier particles produce trails more easily photographed than 
lighter particles. Unless properly corrected, these photographic biases strongly affect 
counts of photographic meteors of any particular type or class. 

For the reasons described, a weighting factor was developed in reference 1 for cor- 
rection of both the velocity and the mass bias. In references 2 and 3 that weighting 
factor was confirmed and refined. A corrected histogram was obtained for velocities of 
meteoroids relative to a hypothetical gravity-free Earth on the basis of data for more 
than 2000 sporadic meteors photographed and reduced by the Smithsonian Astrophysical 
Observatory (ref. 4). A bimodal log-normal equation was derived to express the ve- 
locity distribution shown by the histogram. An older equation for influx rates of parti- 
cles of various masses was confirmed on a higher confidence level than had been possible 
before. 

The velocity histogram and consequently also the bimodal log-normal equation were 
not corrected for spacewise bias. Such a bias exists because particles approaching 
Earth from some directions in space are more likely to enter the atmosphere over New 
Mexico (the camera site) than particles approaching from other directions. Correction 
of spacewise bias would have involved unnecessary complication and the effect of the cor- 
rection would not even have been desirable relative to the primary purposes of the anal- 
yses. In the analysis reported here, use was made of factors for correction of space- 
wise bias for each sporadic meteor, which were calculated as part of the work recently 
reported in reference 5. 

In addition to the more refined results for a gravity-free Earth, similar results 
were obtained for velocity distribution and influx rate relative to the atmosphere of the 
real Earth. The method that was developed for utilization of the results for meteoroid 
damage prediction was also applied to the real Earth condition. 


METHOD OF PROCEDURE 

This analysis covered the following four broad aspects: 

(1) A repetition of part of the effort reported in reference 3, but with a weighting 
factor correcting for spacewise bias in the photography of meteors. This weighting fac- 
tor was in addition to the correction that was made in reference 3 for meteoroid velocity, 
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mass, and other minor parameters. The results were a fully corrected histogram and 
bimodal log-normal equation for particle velocities relative to a hypothetical gravity- 
free Earth (approximately the same as for the real Earth at a distance equal to the ra- 
dius of Earth’s sphere of influence). 

(2) A correction for spacewise bias on the constant factor a, which was determined 
in reference 3, for use in the equation 


F > = am~P (1) 

which, with (3 approximately equal to 1. 34, expresses the influx rate F^ for meteor- 
oids of mass greater than m in grams. (Symbols are defined in appendix A. ) 

(3) Derivation of a histogram and a bimodal log-normal equation for velocities of 
meteoroids relative to the atmosphere of the real Earth. This result was obtained by 
theoretical modification of the histogram for gravity-free Earth conditions, which was 
obtained in the first aspect of this analysis. 

(4) Integration of the bimodal log-normal equations obtained in the first and third as- 
pects of the analysis to obtain average values of velocity to an n 11 power, and an illus- 
tration of how the average value of the n Ln power may be used to calculate necessary 
armor thickness for a space mission. 


Velocity Distribution for Gravity-Free Earth 


The points plotted in figure 1 represent the distribution of Vq, which is the velocity 
of sporadic meteoroids relative to a gravity-free Earth, or the impact velocity on a 
space vehicle moving in Earth orbit but at a great distance from Earth. The data for the 
figure were obtained with the same computer program as described in reference 2 and as 
used in both references 2 and 3. Essentially, that computer program simply counts me- 
teors belonging within each of a large number of intervals of velocity, but with applica- 
tion of a weighting factor to the count of each meteor; that is, if the weighting factor for 
a given meteor has a value <^ w , then such a meteor is counted not as one meteor but as 
<P W meteors. For the data of figure 1, the weighting factor was more refined than as 
used in reference 2 or 3. It is defined by the following equation 


<p w = (cos Z R ) 


- 0.22 


F(Zr) 


0. 712 -3.87 

’oo 

av 


V. 


( 2 ) 


in which Z R is zenith angle, or angle of the meteor path to the zenith, F(Z n ) is a 

K av 

statistical function of zenith angle, position of the meteor within the field of view of the 
cameras, and azimuth of the meteor path, v^ is the velocity of the meteor particle rel- 
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ative to Earth’s atmosphere just before the first deceleration of the particle by atmos- 
pheric resistance, and <p „ is a spatial weighting factor the inverse value of which was 

o 

computed in the effort described in reference 5. The function F (Z R ) was computed 

" n y 

for each meteor as reported and described in reference 1. 

The factor cp differs from the most refined form that was developed in reference 3 
only in the substitution of cp for an approximate correction of the gravity focusing ef- 
feet. A separate adjustment for the effect of gravity focusing was not used (as in refs. 2 
and 3), because the factor cp included an accurate adjustment for that effect. The as- 

O 

sumptions made and the principles followed in calculating the inverse value of cp for 
each meteor as part of the effort reported in reference 5 are described in appendix B. 

From more than 2000 sporadic meteors reported in reference 4, only 1282 were 
chosen for use in obtaining the data of figure 1. The meteors selected for use all had 
values of cos Zj^ at least as high as 0. 2, to avoid unrealistically high values of <p^. 
They all arrived at Earth’s surface from directions in space lying within quadrant 1 as 
illustrated in appendix B. As discussed in that appendix, the velocity distribution for the 
meteors from quadrant 1 may reasonably be assumed to be the same as for all quadrants 
until contrary evidence is available. Moreover, the velocity distribution for quadrant 1 
should more nearly represent the distribution for all quadrants than would a distribution 
obtained with indiscriminate use of all sporadic meteors reported in reference 4. This 
opinion is held because the camera sites in New Mexico could not completely cover any 
quadrant except quadrant 1. 

The curve in figure 1 was fitted to the plotted data. It represents a bimodal log- 
normal equation of the form 


f(v) = C 1 exp<-- 


log e (v + 6) - n 1 


>+ Cg exp 


iAog e v - m 2 n 


(3) 


with values of the constants as follows: 

C 1 = 0.07594 At* =2. 30 <^ = 0.47 5 = 0" 


C 2 = 0.0002317 p 2 =4. 10 ct 2 =0.097 


(4) 


(Equation form (3) contains the constant 6 because that constant will have a nonzero val- 
ue in later use for another purpose. A more usual form of, for example, the first mode 
of eq. (3) would substitute Cj/(v + 6) for C^. The numerical values of and p^ 
would be changed and the physical significance of p^ would be different, but the same 
offset log-normal distribution would be represented. The form used in eq. (3) is more 
convenient for the present purposes. ) In this case, f(v) of equation (3) is f(vQ). The 
function f(v) should be read as ’’frequency of v. ” It may be defined by the equation 
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(5) 


f(v) = 


,[(v 


i dv ) < < ( Tt i dT ), 

dv 


in which p |^v - -i dvj < v r < (v + | dv)] is the probability that an unknown real value 
v r in a specific case lies within a randomly assumed differential range from ~ dvj 
(v + A dv). 


to (v + -t- dv ). Thus, approximately, if Av = 1 kilometer per second is substituted for 

dv and if v is given successive integral values (in km/sec), equation (3) for each suc- 
cessive integral value will give the fraction of all meteors with velocities within 1/2 kil- 
ometer per second of that integral value. 


Revision of Total Unaccelerated Influx Rate 

In reference 3 a value of 2. 98x10"^ g^m'^sec"* was determined for a of equa- 
tion (1). In appendix C, that value is increased to 3. 287x10" ^ because of a correction 
for spatial bias. In consequence, with /3 = 1. 34 as shown in reference 3, it appears that 
the most accurate expression for cumulative mass rate of influx determinable under the 
assumptions and by the methods described in the appendix is 

F > = 3.287X10" 15 m _1 - 34 (6) 

relative to a hypothetical gravity-free Earth. 

Velocity Distribution and Influx Rate Near Earth 

The plotted points in figure 2 represent the distribution of v^ , the impact velocity 
that would exist for a space vehicle in orbit about the Earth at an altitude just sufficient 
to avoid appreciable atmospheric resistance. The data for the plotted points were de- 
rived from those of figure 1 in the manner described in appendix D. The abscissas of fig- 
ure 1 were adjusted for acceleration by Earth's gravity and the ordinates were adjusted 
for the effect of gravitational focusing. Also, increases of ordinates by gravitational 
focusing were integrated to provide a ratio F of influx rate near Earth to influx rate 
far from Earth. The effect of the orbital velocity about Earth for the space vehicle was 
neglected. The statistical effect of the orbital velocity should be small because it would 
subtract from the impact velocity nearly as often as it would add to it and also because it 
would cause a large percentage difference of impact velocity only at the lower impact ve- 
locities, which are not those that represent the principal hazard. In general, the orbital 
velocity of the space vehicle should cause some reduction of the levels of the two peaks 
of the distribution shown in figure 2, some increase in the level of the valley, a decrease 



in the minimum and an increase in the maximum abscissa for which a frequency of ve- 
locity is shown, and some increase in average velocity. 

The curve in figure 2 was fitted to the trend of the plotted points by the same method 
as with figure 1. It represents equation (3) with values of the constants as follows: 

C 1 = 0.1308 p 1 = 1 . 40 0' 1 — 0. 62 5 =-8.0^ 

\ (?) 

C 2 = 8. 966xl0~ 5 p 2 = 4. 14 a 2 = 0. 055 J 

Calculation of Necessary Armor Thickness to Avoid Critical Damage 

Equations (1) and (3), with values of constants given in equations (4) or (7), may be 
used in conjunction with an equation expressing the extent of damage by a meteoroid im- 
pact to determine necessary armor thickness to avoid an unacceptable probability of 
critical damage. The general problem will be analyzed; then a sample calculation will 
be performed for a particular case to illustrate use of the method. 

General expression for critical da mage . - A general expression for a critical value 
of a damage criterion may be derived as follows (under the condition that the criterion 
consist essentially of a constant multiplied by powers of mass m and impact velocity v): 

Let F^r, be the number of impacts per unit area per unit time by particles hav- 
^ ^cr 

ing a combination of velocity and mass such as would cause a damage criterion Z to ex- 
ceed a critical value Z cr - Examples of the damage criterion Z might be necessary 
armor thickness to avoid penetration, or depth of crater formed in a very thick plate. 

In the case of thickness of armor to avoid penetration, the value of Z may be 
(ref. 6) 


Z = c 1 m 1 ^ 3 v (8) 

and, in the case of depth of crater formed in a thick plate, the following equation may 
apply (ref. 6). 

Z - c 2 m 1 / 3 v 2 / 3 (9) 

For the general case, under the condition to which this derivation is limited, 

Z = Cgm^v 6 (10) 

where both A and e are positive constants. In equations (8) to (10), the symbols c^, 


6 
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Cg, and Cg may be constants or functions of variable parameters other than m and v. 

Now, for the general case (eq. (10)) the problem is to consider all meteoroids im- 
pinging on the vulnerable area and to determine which of those meteoroids will produce a 
value of Z greater than Z cr . Actually, no equation is yet available to determine the 
total impingement. However, equation (1) will provide the total impingement for mete- 
oroid masses above a value m x , which may be designated as the minimum value of m 
for which equation (1) applies. For such total impingement, above the mass m , equa- 
tion (1) indicates a value 


F 


>m x 


o»m 


-0 

x 


For any value m above m , equation (1) is equivalent to 

X X 


( 11 ) 


F 


>rn 



dm 


( 12 ) 


where, for m of equation (1), a value m r is substituted. 

Now, from equations (11) and (12), the frequency of m in terms of the total im- 
pingement of meteoroids with mass above m must be 

f(m) = a §_ m ~ (fi +1 ) (13) 

Fv 

>m x 

because then the fraction of the total population under consideration that have mass above 
any given value m r will be 


F 


>m_ 


F 


>m x 



f(m)dm 


(14) 


(Here f(m) is the frequency of m, which may be defined by equation (5) on substitution 
of m for v and m r for v r . ) 

From the definitions of f(v) and f(m), the product f(v)f(m) dm dv must be the frac- 
tion of the total population under consideration that have the mass m ± 1/2 dm and the 
velocity v ± 1/2 dv. So the fraction of the total population that will produce a value of 
Z greater than Z cr must be 


>Z 


cr 


>m. 



f(v)f(m) dm dv 


(15) 
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where m cr , from equation (10), is 


m 


cr 



1/A 


So, from equations (13), (15), and (16) 


>m x >m x y ° ( Z cr/ C 3 v ) 


(16) 


(17) 


with the reservation that the lower limit of the integration relative to m must never be 
less than m x for any practical value of v, that is, for any value of v for which f(v) 
is significantly different from zero. So 


F 


>Z 


cr 


ac § A z-j A r 

y o 


(18) 


or, as the velocity v is always measured positively in the direction of meteoroid motion 
and negative values of v are therefore not real, 


F> z = oc^ X Z~^ X v e/3//A (19) 

(See ref. 8 or other textbook on mathematical statistics. ) Equation (19) may now be con- 
verted to 


Z 


cr 



( 20 ) 


Now, if exposure of a vulnerable area A is required for a time interval T, with an 
expected mean value p for number of impacts involving a value of the criterion Z 
greater than Z cr , then the following equation must apply, 


F 


>Z 


cr 


P 

TA 


(21) 


From equations (20) and (21), with F^, z expressed in the same units as a, 

' cr 

\A//3 


otTAv 


e/3/A 


cr 


( 22 ) 


8 



Equation (22) permits direct determination of Z cr , such as necessary thickness of 
armor, to assure a desired mean number of destructive impacts per mission. 

Negligible error results from using a value of p in equation (22) equal to the com- 
plement of the desired probability p(0) of no destructive impact in a mission, if that de- 
sired probability is 0. 99 or greater. For substantially smaller probabilities of no de- 
structive impact, the following equation based on the Poisson distribution should be 
used (ref. 8 or other textbook on mathematical statistics): 

p = -log e [p (0)] (23) 

Effect of change of influx rate on ex pression for critical damage . - Equation (22) in- 
cludes as factors a value a concerning meteoroid influx rate i n Eart h’s orbit but far 
from Earth and an average value of a power of impact velocity v e £A. In appendix D, 
F c of equation (Dll) is shown to be a factor that would approximately convert the influx 
rate of meteoroids far from Earth to the influx rate near Earth. 

Now equation (22) can be applied to the near-Earth condition by substituting for a 
the product of a and the factor F c and using the average value of the power of impact 
velocity that exists near Earth. However, it will be more convenient for use with the 
velocity distribution near Earth (as shown in fig. 2) and with other velocity distributions 
that might later be found to modify equation (22) to a form in which the constant ol may 
always be used. This objective may be accomplished by combining the factor F c (or its 
equivalent in later cases) with the average value of the power of impact velocity. For 
generality, a function <p„ will be used instead of F which applies specifically to the 
near-Earth condition without regard to direction of impact. The product of a and <p 

O 

may be substituted for a in equation (22) with the result 



For general use, the function ^gV n may be plotted over a wide range of values of n. 

The parameter v in such function could be any one of an infinite variety of impact ve- 
locities. For example, instead of the full velocity on a randomly oriented plane, it could 
be the normal component of impact velocity on such a plane, the normal component of 
impact velocity on a plane with a particular orientation, the normal component on a cyl- 
inder with a particular orientation, and so on. In use of the plot, the exponent n in the 
function <?gV n would be given the value e/3/X as in equation (24), the value of the func- 
tion would be read from the plot for that value of n, and the result would be substituted 
into equation (24). 

For impact velocities without regard to direction of impact, far from Earth, <p is 

O 

equal to unity, as in equation (22), and the impact velocity v of equation (24) is the ve- 
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locity Vq. For the near-Earth condition, without regard to direction of impact, is 
equal to the factor F c of equation (Dll) and the impact velocity v of equation (24) is the 
velocity v^ . 

Sample calculation of critical damage criterion. - A sample calculation will be de- 
scribed here for the purpose of illustration of the method. It is not intended as a defin- 
itive indication of necessary armor thickness, particularly because it involves no con- 
sideration of shower meteors. Because the calculation is intended only to illustrate the 
method, no support will be offered for the values of constants used. For this example, 
the necessary thickness of armor will be calculated for a randomly oriented space radi- 
ator of 2000 square feet vulnerable area, with an assumption of a 99 percent probability 
of no destructive impact in a period of 1000 days. 

Use will be made of an equation by Fish and Summers of Ames Research Center 
(personal communi cation): 


1 = 0. 65 
d 




, 7/8 

norm 


d l/18 


(25) 


where t is minimum armor thickness (cm) at which penetration would occur, d is 
diameter of particle (cm), e. is percent elongation of armor material, p + is density of 
armor material (g/cnr), p p is density of particle (g/cm' 3 ), and v norm is the compo- 
nent of impact velocity normal to the surface of the armor at the point of impact 
(km/sec). With a simplifying approximation that all particles are spherical, equa- 
tion (25) is equivalent to 


t = 0. 8162 e 


-1/18.8/54 -1/2 19/54 0. 875 
t p p p t m norm 


(26) 


(A procedure analogous to the following would apply with any damage criterion that 

might be chosen in preference to equations (25) or (26) so long as a power of v norm 

and a power of m appeared as a factor within that criterion. ) 

Stainless-steel armor will be assumed with values of p t equal to 8 grams per 

cubic centimeter and e^. equal to 10 percent. Meteoroid density p p will be taken as 

0.2 gram per cubic centimeter. On substitution of these values into equation (26), and 

combining of constants, that equation becomes equivalent to equation (10) with use of 

v for v and t for Z, with 

norm ’ 

c 3 = 0. 2000 ^ 

X- 0.3519 f (27) 

e = 0. 875 
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Now equation (24) has been shown to be an expression for a critical value Z cr of 
the damage criterion Z (t in this case) if that damage criterion is expressed by equa- 
tion (10). For the present case, for solution of equation (24), use will be made of a 
value of /3 equal to 1. 34, as was well confirmed in reference 3. Values of c^, A, and 
e from equations (27) will be used. Use will be made of a value of a equal to of 

equation (C8), a value of T equal to the time of the mission (in sec), a value of A 

o 

equal to the vulnerable area (in m ), and a value of p approximately equal to the com- 
plement of the desired probability of no destructive impact 0. 01. The velocity v in 
equation (24) will be taken as v norm , as required for application of equation (26). The 
result for both near-Earth and far -from -Earth conditions and for a critical thickness 


t 


cr 


0. 05046 



3.332\ 02626 
v norm / 


(28) 


In appendix E, equation (Ell) is developed to express the value of v ” orm in terms 
of v n , for any unshadowed randomly oriented surface. From equations (28) and (Ell), 


t 


cr 


0.039005 



^3.332 J 


0.2626 


(29) 


For far -from -Earth conditions in equation (29) cp equals one and v equals v^. For 

near-Earth conditions cp equals F of equation (Dll) and v equals v . 

g c 

For both conditions, values of v n have been determined in manners described in 
appendix F for values of n from 0 to 6, and the results are plotted in figure 3. Also 
plotted in figure 3 are values of F v^, or, for this case, <p v^ . In the calculations of 

o 

the plotted values, equation (3) was used with values of the constants as shown in equa- 
tions (4) and (7). 

3 332 

Values of the function cp_v for far- from- Earth and near-Earth conditions, 

somewhat more accurate than can be read from figure 3, are 


Vq 332 = 18 270 


•h 


2 . 447 v 3 • 332 = 35 030 


(30) 


From equations (29) and (30), for far- from- Earth conditions 


t cr = 0. 513 cm (0. 202 in. ) 


(31) 


and for near-Earth conditions 
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(32) 


t cr = 0. 609 cm (0. 240 in. ) 


In the transformation from equation (17) to (18) a reservation was necessary that the 
expression for m in equation (16) must always exceed the value of m the lower limit 
of mass for which equation (1) applies. From figures 1 and 2, the highest practical value 
of impact velocity is certainly lower than 80 kilometers per second. The lowest prac- 
tical value of m „ will therefore be not lower than 
cr 


m 


cr(min) 


0 513 \ 1/0.3519 _a 

u-0i,i ' = 2. 7x10 g 


^0. 2x80 

for the condition far from Earth, and 


0. 875 


m 


cr (min) 


^0. 2x80 


0 609 \l/0. 3519 a 

) =4. 4X10 4 g 

0. 875 / 


(33) 


(34) 


for the condition near Earth. By reference to figure 7 of reference 3 and the pertinent 

-4 

textual discussion, m appears to be at least as low as 6. 3x10 . Moreover, a 

x -4 

straight-line extrapolation to the value 2. 7x10 appears to be quite reasonable. Hence, 
the reservation appears to be satisfied. 

For comparison, the same problem will now be solved by a method heretofore in 
common use with an average impact velocity of 30 kilometers per second that has also 
been in common use for substitution into equation (10). Then, with the values shown in 
equation (27) and with substitution of a critical thickness t for Z, equation (10) yields 
a value of critical mass 

. 1/0.3519 

m cr = ( = 0. 02058 t^ 841 (35) 

'0. 2x3 0 0, 875 / 


Next, with substitution of m cr for m and p for F^, in equation (6) and multiplication 
of the result by TA, 


p = 3.287X10" 15 x (0.02058 t^ 841 ) 8. 64X10 7 x 185.8 

= 9. 6016X10" 3 t -3- 807 


(36) 


Finally, for p = 0. 01, 


t cr = 0. 9894 cm (0. 3895 in. ) 


(37) 
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The critical thickness of 0.3895 is nearly twice the value shown by equation (31) for 
the unaccelerated condition. However, if the older method of calculation in common use 
is given the benefit of using only the normal components of velocity, the contrast of re- 
sults is far less striking. Thus, by equation (Ell), the assumed average value v^ = 30 
yields an average of the normal component VQ^ norm ^ = 20. Use of this value in equa- 
tion (35) produces a value of t equal to 0. 2732 instead of the value shown in equation 
(37). This value is only 35 percent greater than the value shown in equation (31). 

A much greater reduction of critical thickness might superficially be expected be- 
cause the average impact velocity is now known to be less than one half the value of 
30 kilometers per second previously assumed and used. The much greater reduction did 
not materialize because the older method erred not only in the use of an average velocity 
much too high, but also in its assumption that the critical thickness would even depend on 
the simple average velocity. As may be seen from equation (29), for the damage crite- 

rion used, the critical thickness really depends on the function (v 0- 4 ) . Unfor- 

tunately, the function v n grows more rapidly than does (v) n as n increases above 

unity. From figure 3, the value of (v J ‘ J ^) is approximately 12. 9. By the older 

method, with use of equation (26), the corresponding value would be 30^' or 19. 6 . 

3 33 2 0. 2626 

This value is only 52 percent greater than the value of (vq ) . This comparison 

would call for 5 2 perc ent greater critical thickness by the old method. However, with the 
older method, v ^ orm according to equation (Ell) is a smaller fraction of v n than with 
the method presented here, with the result that the older method calls for only the 
35 percent greater thickness noted. 


DISCUSSION OF RESULTS 

Equation (3), with values of constants as in equations (4), appears to be the most 
accurate that can be derived, by the methods used, for the distribution of velocity Vq. 
Such distribution is without regard to direction of impact; that is, it is the expected dis- 
tribution of impact velocity on a randomly oriented surface. It differs from the final 
equation developed in reference 3 only because of correction for spatial bias. The histo- 
gram is but little different than that obtained in reference 3, but it permitted a more ex- 
act fit of a bimodal log-normal equation. Because the photographs were taken at all fea- 
sible hours of the night and on dates well distributed throughout the year, a degree of 
statistical masking of the effect of the spacewise bias was expected. The slight un- 
masked bias shown by the results may exist for the following reasons: ( 1 ) photographs 
were exposed predominantly during later hours of the night, a condition that favored me- 
teors in retrograde orbits, ( 2 ) the photographic site was north of the equator, a condition 
that favored meteoric particles approaching Earth from directions north of the ecliptic 


13 



plane, and (3) the photographs were necessarily exposed entirely at night, a condition 
that favored meteoric particles approaching Earth from the direction opposite the sun. 

Equation (3), with values of constants as in equations (7), seems to be the most ac- 
curate that can be derived, by the methods used, for the distribution of velocity v TO . It 
also is without regard to direction of impact. 

A value of a for use in equation (1) equal to 3. 287x10” ^ g^m’^sec"'*' {a^ ^ of 
eq. (C8)) appears to be the most accurate determination that can be made by these meth- 
ods. That value applies for statistically expected impacts on a randomly oriented sur- 
face moving in Earth’s orbit, but not near enough to Earth that Earth’s gravitation can 
appreciably increase the frequency of impacts. At Earth's atmosphere, the value of a 
should apparently be increased by the factor F c of equation (Dll). 

The method of calculation of the critical value of a damage criterion, with use of 
armor thickness to avoid penetration as a specific example, is a culmination of all the 
effort reported in references 1 to 3 and in the present analysis. It appears to rest on a 
good foundation, but it is incomplete in that shower -associated meteoroids are neglected. 
Ultimate modification to include the effect of shower-associated meteoroids would need 
to include at least one additive term on the right side of equation (12) and a repetition of 
all work on which equation (12) has an effect. Moreover, in equation (3) and all work af- 
fected by it changes would have to be made to account for the frequency f(v) applying to 
shower-associated meteoroids. 

The additional hazard due to shower-associated meteoroids may be great or small, 
depending on the directionality pattern of sporadic meteoroids. Considering as a simple 
example a single radiator tube, later study may show that a much different hazard due to 
sporadic meteoroids exists when the axis of that tube points in the direction of the apex 
of Earth's motion than when it points toward the sun, or in a direction normal to the eclip- 
tic plane. In that event, the optimum result might involve one of two recourses: (1) use 
of the preferred alinement with additional thickening of tube walls to protect against the 
additional hazard from shower-associated meteoroids and (2) alinement of the tube axis 
parallel to the paths of shower-associated meteoroids while passing through showers, 
with additional thickening of tube walls to protect against the additional hazard from spo- 
radic meteoroids while the tube axis is so alined. 

On the other hand, if later study should show little effect of directionality on the 
hazard from sporadics, then the axis of the tube would always be alined parallel to the 
paths of shower-associated meteors. The only additional hazard to the cylindrical part 
of the tube wall would exist if the tube passed through two showers simultaneously. 

Even aside from any uncertainty regarding the effects of meteor showers, no argu- 
ment will be presented here for the exactness of the results as shown by equations (31) 
and (32). The calculation is an illustrative example only. If the reader prefers any 
other criterion than that shown by equation (25) or if he prefers other values of constants 
than those used, he may quickly arrive, by the same methods used herein, at values 
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other than those shown by equations (31) and (32). A part of the purpose of this paper is 
to enable him to do so. 

The more exact results that have been developed in references 1 to 3, and in the 
work reported here, have greatly reduced the average velocity of meteoroids below ear- 
lier estimates. The consequent decrease in necessary armor thickness is not great be- 
cause of compensating inaccuracies in the earlier treatment. But the new, slightly 
lower, calculated thickness is believed to be far more firmly grounded. The compensat- 
ing inaccuracies might have proved serious if a much lower average velocity had not 
been found. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, July 14, 1969, 

120-27. 
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APPENDIX A 


SYMBOLS 


A vulnerable area of armor, nr 

C with subscripts 1 or 2, constant factor for mode 1 or 2 of bimodal log- 

normal velocity distribution 

c unsubscripted, constant factor for a unimodal log-normal velocity distribu- 

tion or for unspecified mode of a bimodal log-normal distribution; with 
subscript 0, 1, 2, or 3, factor containing constants and variable parame- 
ters other than meteoroid mass and impact velocity within various damage 
criterions Z (eg for general case) 

d diameter of meteoroid, cm 

F value of cp for influx at Earth's atmosphere without regard to direction 

c g 

F(Z d )„ average statistical function of zenith angle Z,-,, position of meteor within 

±1 aV Xv 

camera field of view, and azimuth of meteor path (ref. 1) 

F^ number of meteoroid impacts from all directions of mass greater than a 

given lower limit at distance from Earth equal to Earth's sphere of influ- 

- 2-1 

ence, m sec 

Fv^ number of impacts for which damage criterion Z will exceed the value 

Z cr , m^sec- 1 

f( ) with variable argument within parentheses, frequency of a given value of 

argument, defined for argument v by equation (5) 

m mass of meteoroid before any ablation by atmosphere, g 

m cr meteoroid mass at a given impact velocity just sufficient to cause damage 

criterion Z to exceed the critical value % cr , g 

p probability of destructive impact during a mission 

p( ) with argument within parentheses, probability of the argument 

T duration of mission, sec 

t maximum thickness of armor at which penetration would occur on impact by 

a given meteoroid, cm 

tl,t2,t3 (subscripts) first tentative, second tentative, third tentative 

t critical thickness of armor that will provide a specified probability of pene- 

tration by a meteoroid, cm 
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V G 

v norm 

v„, 

v oo() 

v -(v G ) 


Z 


Z cr 

Z R 


B 

5 


c 

e t 

A 


meteoroid velocity oi nature to be specified wherever used, km sec - A 

velocity of escape from Earth's gravitational field at position 90 km above 
Earth's surface, km sec”* 

velocity of meteoroid relative to Earth at distance equal to radius of Earth's 
sphere of influence, km sec”* 

component of impact velocity normal to surface at point of impact, km sec”* 

velocity of meteoroid relative to Earth just before deceleration by Earth's 
atmosphere begins, km sec”* 


with number as argument within parentheses, same as 
equal to argument, km sec"* 



with 


v 


G 


velocity of meteoroid relative to Earth at entry into Earth's atmosphere 
after acceleration from velocity v^ relative to Earth at a distance equal 
to radius of Earth's sphere of influence, km sec”* 


damage criterion relative to meteor impact, such as minimum thickness of 
armor to avoid penetration or depth of crater formed 

critical value of damage criterion Z at which probability of destructive 
damage is just tolerable 


angle of meteor path to zenith, or zenith angle 

constant factor governing rate of meteoroid impact in eq. (1), g H m sec 

same as a but restricted to meteoroids arriving from quadrant 1 as de- 
ft _2 - 1 

scribed in appendix B, g H m sec 

flux rate for impacts of meteoroids on unit surface that is arbitrarily 

turned to orientation perpendicular to path of each approaching meteoroid, 
unspecified units 

expected flux rate for impacts of meteoroids on randomly oriented unit 
surface, unspecified units 

negative value of exponent of mass in equation for cumulative influx rate 
relative to mass 


in log-normal distribution equation for function f(x), offset from argu- 
ment x to be used in measuring distances along logarithmic scale 
from modal position 

exponent of impact velocity in damage criterion Z 

percent elongation of armor material 

exponent of meteoroid mass in damage criterion Z 
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unsubscripted or with subscript 1 or 2, modal value in unimodal log-normal 
equation or value of mode 1 or 2 of bimodal log-normal equation 

-3 

density of meteoroid, gem 

-3 

density of armor material, g cm 

unsubscripted or with subscripts 1 or 2, standard deviation in unimodal 
log-normal equation or in mode 1 or 2 of bimodal log-normal equation 

factor by which a may be multiplied to provide analogous influx rate 
under specified condition other than that to which a applies 

weighting factor for use in count of meteors to correct for biasing effect of 
direction in space from which meteoroid arrives 

weighting factor applied in count of meteors to correct for biasing effects of 
zenith angle Z^, function F(Z^) av , meteoroid mass m, meteoroid ve- 
locity relative to Earth's atmosphere v M , and direction in space from 
which meteoroid arrives 

weighting factor to apply to count of meteors to correct approximately for 
effect of gravity focusing 



APPENDIX B 

WEIGHTING FACTOR FOR CORRECTION OF SPACEWISE 
BIAS IN METEOR PHOTOGRAPHY 

The values of cp as used in equation (2) are based on work by Albers and Diedrich 
(ref. 5). With the exception of meteors having a cosine of zenith angle less than 0.2, 
they computed the inverse value of a weighting factor for each sporadic meteor reported 
in reference 4 to correct for the biasing effect of the direction in space from which such 
meteor came. Certain preliminary assumptions made possible the definition of a set of 
coordinate axes about which the meteoroid flux should be symmetrical. The symmetry 
allows a logical correction for spacewise bias. 


Preliminary Assumptions and Resulting Set of Coordinate Axes 

The following two assumptions were made: (1) For each meteoroid arriving at Earth 
from a direction north of the ecliptic plane, an equal probability exists that a meteoroid 
might arrive in a symmetrically oriented orbit from a direction south of that plane, and 
(2) for each meteoroid arriving at Earth from the spatial hemisphere opposite the sun, 
an equal probability exists that a meteoroid might arrive in a symmetrically oriented 
orbit from the spatial hemisphere toward the sun. These assumptions give full scope to 
(1) the tendency of meteor orbits to lie within planes that make only small dihedral an- 
gles with the ecliptic plane, (2) the tendency for some values of eccentricity of meteoroid 
orbits to be more frequent than others, and (3) the tendency for some values of aphelion 
or perihelion to be more frequent than others. Equivalent assumptions would be that 
(1) for a given inclination of meteoroid orbit to the ecliptic, a meteoroid is as likely to 
intersect Earth's orbit at the descending node as at the ascending node, (2) equal proba- 
bility exists for a meteoroid to pass through Earth's orbit during its transit from aphe- 
lion to perihelion as for it to pass through Earth's orbit during its transit from perihe- 
lion to aphelion. 

The two assumptions may be discussed with reference to sketch (a), which illus- 
trates a coordinate system based on the symmetry involved in those assumptions. The 
origin of the set of axes x-x, y-y, and z-z, as seen in the sketch, should be understood 
as always moving along with the Earth, and the x-x and y-y axes should be understood 
as always rotating about the origin in such manner that the positive direction of the x-x 
axis always points in the direction of the apex of the Earth movement and the positive di- 
rection of the y-y axis always extends in the direction within the plane of the ecliptic, 
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(a) 

perpendicular to axis x-x, and approximately in the direction toward the sun. 

All four quadrants illustrated in the sketch are associated with the two assumptions 
mentioned. Those assumptions, with reference to the sketch, are equivalent to the 
assumption of two planes of symmetry: (1) the ecliptic plane, that is, the plane of 
the x-x and y-y axes, and (2) the plane tangent to Earth’s orbit and perpendicular 
to the ecliptic plane, that is, the plane of axes x-x and z-z. More specifically, 
for every meteoric particle intersecting Earth’s orbit at any angle a to the x-z plane 
and within quadrants 1 or 4, an equal probability exists for a particle from quadrants 2 
or 3, with path at the same angle a to the x-z plane and with the same projection of 
path on the x-z plane. Also, for every particle intersecting Earth's orbit at any angle 
6 to the x-y plane and within quadrant 1 or 2, an equal probability exists for a particle 
from quadrant 4 or 3, with path at the same angle 5 to the x-y plane, and with the 
same projection of path on the x-y plane. 

The advantages of the use of the coordinate axes shown in the sketch, with the two 
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planes of symmetry, will now be discussed. The degree of reasonableness of the as- 
sumptions will be considered somewhat later. 


Advantages of Use of Four Symmetrical Quadrants of Space 

A meteoric particle arriving at the Earth orbit from any direction whatever within 
quadrant 1 could possibly be photographed by the cameras in New Mexico at least at 
some hour during the total time throughout which they were operated in obtaining the pho- 
photographs reported in reference 4. But particles arriving from any one of many di- 
rections within each of quadrants 2, 3, and 4 would have no possibility of being photo- 
graphed at any time. With a suitable spatial weighting factor, and with a large enough 
sample, determination of a complete directional distribution would be possible for quad- 
rant 1, though not for quadrant 2, 3, or 4. A directional distribution was not sought in 
this analysis. However, the velocity distribution that was sought is so closely associ- 
ated with direction of meteoroid travel that complete directional coverage is necessary 
in order to allow a complete determination of the velocity distribution. 

One of the advantages of dividing the space surrounding Earth into the four quad- 
rants as shown in the sketch is now apparent. Meteors arriving from quadrants 2, 3, 
and 4 may be eliminated from consideration for statistical purposes. Then, with a 
proper spatial weighting factor for quadrant 1, an accurate velocity distribution may be 
obtained for that quadrant including particles arriving from any possible direction within 
the quadrant. Because of symmetry, the result should also apply to quadrants 2, 3, and 
4, and, hence, for all space surrounding the Earth. 

An additional advantage is concerned with determination of the total influx rate of 
meteoroids, that is, the value of a in equation (1). The site of the cameras in New 
Mexico cannot be assumed to encounter the same influx rate of particles as the statistical- 
ly expected rate for a randomly oriented surface. The fact that the cameras were oper- 
ated mostly during the later hours of the night favors impact by the comparatively rare 
particles in retrograde orbits, but probably disfavors impact by the average particle in a 
direct orbit. Location of the cameras nearer the equator than the pole tends to favor im- 
pacts by particles in orbits having only small inclinations to the ecliptic plane. It is well 
known that orbits having small inclinations to the ecliptic are more numerous than those 
having larger inclinations. These uncertainties can be eliminated by restricting the 
analysis to meteors arriving from quadrant 1 and with use of a spatial weighting factor 
that relates the probability of an impact on the atmosphere over New Mexico to the prob- 
ability of impact on an equivalent area randomly oriented. The properly corrected influx 
rate for quadrant 1 alone needs only be multiplied by four to obtain a corrected influx 
rate for all quadrants. 
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Justification of Assumptions 


The assumptions of symmetry that have been used are a necessary corollary of a 
concept of randomness of orientation of meteoroid orbital planes, namely, that the ori- 
entation is subject to a distribution law favoring small inclinations to the ecliptic plane 
but that the orientation is otherwise random. Such a condition is not known to exist but 
has also not been disproved. Because of the impossibility of photographing meteors in 
daylight, symmetry between quadrants 1 and 2 or between quadrants 3 and 4 could be 
checked eventually only by methods other than photographic. Symmetry between quad- 
rants 1 and 4 could be checked by repetition of the work of references 1 to 4 and by rep- 
etition of the analysis reported here, with use of cameras located within the southern 
hemisphere. Until such checks become possible, the assumptions of symmetry appear 
to represent the simpler condition and should be used until discredited. 

A definite possibility exists that meteoroids not identified with known showers are 
all nevertheless members of unknown minor showers. These minor showers might pos- 
sibly be roughly grouped in a manner such that the sporadic meteors could be regarded 
as members of a few very diffuse showers. "Diffuse shower, " here, is a shower in 
which the true radiants of the different meteoroids are not identical but are more or less 
approximately the same. A possible consequence would be that the sporadic meteors 
would tend to vary in the distribution of their directions from one part of the year to 
another as expressed relative to the coordinate axes that have been illustrated. The 
correct integrated effect throughout the year, then, would differ from that obtained in 
reference 5 because of the tilt of the Earth's axis relative to the ecliptic north. It is not 
believed that the error would be great. 


Spatial Weighting Factor 

In addition to the determination and recording of the inverse value of a spatial 
weighting factor for sporadic meteors reported in reference 4, Albers and Diedrich 
(ref. 5) determined and recorded the quadrant from which each meteor arrived. For 
their use, the authors of reference 4 supplied data as to the number of minutes of each 
hour of each day that the cameras were operated throughout the period covered by the 
data reported in reference 4. 

For each meteor, Albers and Diedrich performed the following operations: (1) For 
each hour of camera operation, a ratio was determined equal to a unit cross-sectional 
area perpendicular to the meteor path at a great distance from Earth, through which the 
meteor might have passed, divided by the corresponding area on the surface of the Earth 
at the camera sites in New Mexico within which the meteor might have arrived. The de- 
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termination of this ratio included a determination of the hyperbolic orbit of the particle 
about Earth center and an accurate determination of the effect of gravitational focusing. 
(2) The ratio so determined for each hour was multiplied by the number of minutes of 
camera operation within that hour. (3) All results of item (2) were added and divided by 
the total number of minutes of operation of the cameras. The reciprocal of the result is 
the spatial weighting factor cp s for the meteor as used in equation (2). 
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APPENDIX C 


REVISION OF INFLUX RATE FOR GRAVITY-FREE EARTH 
FOR CORRECTION OF SPATIAL BIAS 

The value of 2. 98x10“^ g^nT^sec"* that was found in reference 3 for a of equa- 
tion (1) will be modified to correct for the spacewise bias caused by the fact that an area 
in the sky above the camera location would not necessarily receive the statistically ex- 
pected frequency of impacts for a similar area located at random over Earth's surface. 

- 15 

Determination of the factor 2. 98x10 involved use of a basic reference value from 
reference 9 and an actual count of sporadic meteors of mass greater than 0. 705 gram. 
The count was corrected for an estimated fractional failure of observation because of 
faintness of meteor traces and approximately corrected for gravitational focusing by ap- 
plication of the factor 


(C1) 

where v g is the velocity of escape from Earth's gravitation (ref. 10). The lower mass 
limit of 0. 705 gram was an arbitrary value above which observational failure caused by 
meteor faintness was believed to be small. Because it was believed to be small and be- 
cause the failure should apply to all spatial quadrants shown in appendix B about equally, 

correction of this observational failure will be ignored in the following treatment. Ac- 

-15 

cordingly, the value of a, 2. 98x10 , may be reduced to a value applying to quadrant 1 

alone by the equation 


K, 


a l ~ 


K 


a 


(C2) 


1-4 


where K-^ is a count of sporadic meteors within quadrant 1 and Kq_ 4 is a count of spo- 
radic meteors for all quadrants, in both cases with cos Z R not less than 0.2, with 
mass greater than 0. 705 gram, and with correction for gravitational focusing by applica- 
tion of the factor <P w (g) °f equation (Cl). The corrected counts Kj and Kj ^ and the 
consequent value of the ratio K^/Kj_^ were determined by a minor addition to the com- 
puter program used for the results described in the preceding section. The result was 


K, 


= 0.608 


K 


(C3) 


1-4 
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or 


a 1 = 1. 812x10" 15 (C4) 

Now the adjusted count of equation (C2) may be expressed as 

K 1 * E^w(G) < c5 > 

where the summation is performed for all sporadic meteors within the classification de- 
scribed, with each meteor counted as <P w (g) meteors. This adjusted count may be fur- 
ther adjusted for spacewise bias as 


K 1( S )=I>S 

and an adjusted value of a may be obtained as 


“l(s) 


K 


l(s) 




K, 


(C6) 


(C7) 


a determination of which was included in the computer program, with the result 

Q! 1(s) = 1. 814 oi x = 3.287X10" 15 (C8) 

The value a i( s ) is applicable to sporadic meteors arriving from quadrant 1 only. 

The factor of course, refers to impacts occurring on the atmospheric surface as it 

existed over the camera sites. But, as the factor cp is referred to an area arbitrarily 
alined normal to the path of each meteor, the factor Q^g) also applies under the as- 
sumption that the vulnerable area would always be so oriented as to be normal to the path 
of any approaching particle. As shown in appendix E (eq. (El)), a randomly oriented 
surface should receive only one-quarter as many impacts. However, applies to 

one quadrant only. By the assumption of symmetry, the total flux from all quadrants 
should be four times as great. The two adjustments cancel, so that the value of 
shown in equation (C8) represents the statistically expected total influx from all direc- 
tions on one side only of a randomly oriented surface. Hence, this value is used in equa- 
tion (6) for the cumulative mass distribution. 
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APPENDIX D 


DERIVATION OF VELOCITY HISTOGRAM AND INFLUX 
RATE FOR NEAR-EARTH CONDITIONS 

The histogram represented by the circular points in figure 2 was derived theoreti- 
cally from that of figure 1. Use was made of the relation, for change of variable, 

f(vj = f(v Q ) dv G /d Voo (Dl) 

(see ref. 8 or other textbook on mathematical statistics). For application to a histogram, 
equation (Dl) must be converted to an approximate form 

f(vj = f(v Q ) Av g /Av to (D2) 

Equation (D2) converts from frequency of v G to frequency of v w , but only with regard to 
the relative values of the two variables. Account must still be taken of the effect of 
Earth's gravity in focusing the influx of particles in such manner as to increase the rate 
of influx for all particles, and to increase the rates of influx for slower particles relative 
to those for faster particles. Thus, equation (D2) becomes 

f( v J = C g f ( v G ) Av g /A v °° (D3) 

where is the concentration factor for the given value of v G due to the gravitational 
focusing effect. Finally, the values of f(v 00 ) must be normalized to provide a total value 
of unity. 

Details of the method of deriving the histogram shown in figure 2 in accordance with 
equation (D3) will be described by considering an illustrative example, as follows: 

In the data for figure 1, four values of f(v G ) are 

f(l. 5) - 0.002704' 
f(2. 5) = 0.01030 

r (D4) 

f (3. 5) = 0. 005130 
f (4. 5) = 0. 02960 ^ 

Now, for each of the frequencies of v G shown in equations (D4), an augmented velocity 
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was determined by use of the equation 


( y G ) 


v -(v G )= V v G + v e 

2 

where v g is the square of the escape velocity, taken as 

v 2 = 123.25 km 2 /sec 2 


(D5) 


(D6) 


The left side of equation (D5) may be regarded as an argument of the frequency function 
at the left side of equation (D3). The augmentation of velocity from Vq to ^ with 

use of equation (D5) was based on the obvious fact that the kinetic energy of a particle 
must be increased by the energy of escape when a particle moves from an infinite dis- 
tance to the Earth's surface. Very little error is involved, even at low values of Vq be- 
cause of the incorrect treatment of particles as being accelerated from an infinite dis- 
tance. The results of the augmentation of velocity for the four examples were 


Ml. 5) 

M2.5) 
M3. 5) 
v °°(4. 5) 


11. 203 > 
11. 380 

► 

11.640 
11 . 979 , 


(D7) 


The frequencies of Vq, that is, f(v G ) as given by equations (D4), represent influx 
rates at the values of Vq given (±0. 5 km/sec), divided by the total influx rate far from 
Earth for all velocities. Tentative values of the same frequencies augmented by gravi- 


tational focusing f 
and 10) 



were obtained with use of the following equation (refs. 2 


f 




(D8) 


The left side of equation (D8) corresponds to the left side of equation (D3). The right 
side of equation (D8) corresponds to the first two factors of the right side of equation (D3), 
but not yet in their final form. The results from equation (D8) are normalized relative 
to the total influx far from Earth rather than relative to the total influx near Earth. Also, 



represents a fractional influx rate within an interval of 1 kilometer per sec- 
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ond in the value of v G , not the value of v m 


< V G>' 


Equation (D8) involves substantial error at low values of v G because its derivation 
assumes that the meteoroid is accelerated by Earth gravity from an infinite distance, 
while, in fact, because the particle was originally in orbit about the sun, its acceleration 
by Earth gravity does not occur from an infinite distance. An estimate of this inaccuracy 
was presented in figure 8 of reference 2. From that figure, for each value of v G used 
in equation (D8), a correction factor cp er was read and used as follows to obtain a sec- 
ond tentative frequency value 


ik 


>(v r ) 


= f 


t2 


Kv-) 


<P 


tl 


er 


The results, for the four cases of equations (D4) and (D7) were 

* [mi. «>L-°- 09426 

f [ v ~( 2 . S )] t2 = ° 1905 
f [ v ~(3-5)] t 2 = °- 05328 
£ [ V -( 4 . 5)] t2 =°- 2021 „ 


(D9) 


(DIO) 


The left side of each of the equations (DIO) again corresponds to the left side of equa- 
tion (D3), while the right side of these equations corresponds to the first two factors in 
the right side of equation (D3) in their final form. At this stage, the values of 


3 (v r ) 


were still in units of fraction of total flux at great distance from Earth per 


t2 


unit value of Av, 


The conclusion follows that the total of all values of f l v °°(v G ) 


■ v, 


ob- 


Jt2 


tained with use of equation (D9) is a concentration factor, equal to the ratio of total flux 
rates per unit area at Earth's atmosphere and at a great distance from Earth but within 
the Earth orbit. The value of such concentration factor proved to be 


F c = 2. 447 


(Dll) 


The results of the use of equation (D9), as shown in equations (DIO), were next con- 
verted to fractions per unit value of , % , but still normalized relative to the total 

' G' 

influx rate far from Earth, with the equation 
28 



f |" v “(v^)l =f [ v “(v„)l x “ 

L ' v G)J t3 L ' G'J^ 4 v„ (V q) 

in which 

Av °°( v g ) = 2 [ V °°( V G +1 ) " V °°( V G -1)] 
except in end cases, in which 

Av °°(v g ) = V °° (v G + 1 ) ■ Mv Q f 
or l 

AV °°( v g) =Vo0 ^ " Voo ^ v G -1 ^ 


(D12) 


(D13) 


(D 14) 


Equation (D12) yields a final equivalent result to equation (D3). Equations (D13) and (D14) 
yield approximate averages of the intervals between one value of v 


Hence, f 


d (v g ) 

^ Jt3 


°°( v g) 


and the next. 


given by equation (D12) is a fraction of total influx far from Earth at 


the near Earth velocity ^ ±1/2 kilometer per second. 

The results of application of equations (D12) to (D14) to the values shown by equa- 
tions (D7) and (DIO) are 


Av oo(l. 5) - v °°(2. 5) " Ml. 5) - 11 * 380 " 11 ■ 203 “ 0> 177 
Av »(2.5r}[ T ”(3-5>‘ V “>(l-5>] <11.640 - 11.203) = 0.2185 > 

A '»(3.5) = ~ ['-(4.5) - '-(2.5)] = ~ tt 1 - 979 - U- 380 * * °- 2995 _ 


(D15) 


and 
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1 

_ 0. 1905 = 

0.872 l 

[ V “(2. 5) 

t3 


v °°(2. 5) 

t2 

Av t 

*(2. 5)_ 

0.2185 

[ V -(3. 5) 


- f 




1 

_ 0. 05328 

= 0. 178 

t3 

— 1 

v °o(3. 5) 

t2 

Av t 

»(3.5)_ 

0.2995 

J 


(D16) 


The frequency values shown by equations (D16) and the other similar frequency val- 
ues not shown were still not normalized relative to the total near -Earth flux. 

Finally, for values of v^ equal to 1. 5 , 2 . 5 ,... 79. 5 , frequencies fCv^) were ob- 


tained by linear interpolation relative to the values of v 0 


f 




>(v„) 


and the values of 


of equations (D7) and (D16). The results for values of v M from 1. 5 to 10. 5 


t3 


were taken as zero because these values of 
interpolation was executed between the values of v 
relative to the corresponding values of f Jv^g 5)J 

The unnormalized result of this interpolation is 


v^ were less than 


°°(2. 5) 
and f 
t3 


and v. 


3 (1- 5)' 


For v = 11.5, 


(3.5) 


_ V °°(3. 5)] 


t3 


of equations (D7), 
of equations (D16). 


f (11- 5 ) un = °- 552 (D17) 

The result as shown by equation (D17) and the similar results for all values of 
from 12. 5 to 79. 5 were normalized by dividing each value by the total of all, which was 
2. 571. Thus, for the case of v^ equal to 11. 5 kilometers per second, the normalized 
result was 


« U - 5 >norm" 0 - 2147 < D18 > 

The normalized results for all velocities from 11.5 to 79. 5 kilometers per second, con- 
verted to percentages, are plotted as the circular points in figure 2. 
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APPENDIX E 


AVERAGE n th POWER OF NORMAL COMPONENT 
OF IMPACT VELOCITY 

As part of the problem of necessary armor thickness for meteoroid protection, de- 
tailed consideration must be given to the statistical relation between normal components 
of impact velocity and the full impact velocity for a plane surface exposed on one side to 
an isotropic flux of meteoroids. The problem is identical to that involving the statisti- 
cally expected flux arriving from one direction only, or from any distribution of direc- 
tions, on one side of a randomly oriented surface. In one case, it is the direction of 
movement of particles that varies randomly while the surface remains stationary. In the 
other case, the surface orientation varies randomly while the flux remains fixed. 

This problem as stated is independent of the nature of the velocity distribution, 
whether log-normal or of some entirely different nature. For the case of the plane with 
fixed orientation, the only assumption is that the meteoroid flux is isotropic. For the 
expected relation between normal component of velocity and the full velocity, in the case 
of the randomly oriented plane, even the assumption of isotropicity is not needed. 

At the outset, all meteoroids will be treated as having the same velocity v, relative 
to the point of impact on the s urface . Later, a velocity distribution of unspecified nature 
will be considered. A value v^ orm will be sought, an expected arithmetic average of 
v norm’ w ^ ere v norm * s com P onen t of the impact velocity normal to the surface. 
Impacts of any obliquity whatever will be included in the treatment. 

The following proposition will prove useful: given an influx rate of meteoroids a 
measured as impacts per unit area per unit time on one side of a plane surface that is al- 
ways arbitrarily oriented to face any approaching particle, prove that one side of a ran- 
domly oriented surface will receive an influx at r = (1/4)0?^. 

Now, a , as defined to represent the given condition, is the same as the average 

ir 

number of impacts of any obliquity whatever per unit time on the total surface of a sphere 
per unit area of a cross section constructed through the center of the sphere. This fact 
may be seen by the consideration that a sphere offers exactly its cross sectional area for 
interception of a meteoroid approaching from any direction. But a r as defined (the in- 
flux on one side of a randomly oriented surface) is the same as the average number of 
impacts per unit time on the total surface of a sphere per unit area of that surface. This 
fact is apparent in view of the consideration that a sphere has infinitesimal elements of 
surface facing impartially in all possible directions. Hence, over an infinite span of time, 
the statistical effect of the infinitesimal elements would be the same as that of a plane 
area that changes its orientation randomly from time to time. (Every impact will involve 
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two punctures of a sphere, one on entering and one on leaving. The two punctures, for 
the purpose of this discussion, constitute one impact. ) 

From the foregoing considerations 


a 


r 



(El) 


Now detailed consideration will be given to isotropic impacts on one side of an infini- 
tesimal plane area A as shown in sketch (b): 


z 



For this purpose, a hypothetical hemispherical surface of radius p will be considered, 
with the infinitesimal area A at its center. Equal influx rates should exist through all 
parts of that hemispherical surface, passing through the elementary surface A, if that 
elementary surface is always oriented to face any approaching particle. A differential 
element ds of the hemispherical surface will be considered, bounded by two circles 
centered about the axis z-z and separated by an arc d 6 (where 0 is the angle from the 
z axis). 

Now an expression d cp will be sought for the flux-area product (for the infinitesimal 
area A) passing through the surface ds and impacting on the surface element A. The 
flux-area product for impingement on the area A from any direction whatever, if the 
area A were always oriented normally to the path of any approaching particle, would be 
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(E2) 


d< ?l = a p A = 4a- r A 

Under the assumption of isotropicity, the part of that flux-area product that would pass 
through the annular surface ds would be 

2 

d cp 2 = d <p^ 2 — sin 6 d - ■ = 2a r A sin 6 dd (E3) 

4 up 2 

But the area A is always oriented at an angle 0 to the flux passing through the annular 
surface, so that the flux-area product becomes 

dcp = cos 0 = 2a? r A sin 6 cos 6 d0 (E4) 

The factor A cos 6 in equation (E4) represents the projection of area A on a plane nor- 
mal to the meteoroid path, that is, the effective area presented by A for interception of 
meteoroids approaching at the angle 9. 

in 

Now the arithmetic average of the n n power of normal component of impact velocity 
(v constant) may be expressed as 


v 


n 

norm(v=k) 


rir/2 

Je^o 


v norm(v=k) dcp 


« r A 


(E5) 


The numerator of equation (E5) is the integrated value of v n 0rm ( v= k) multiplied by 
the differential flux-area product. The denominator, by previous definition, is the total 
flux-area product for the one side of the area element A. 

From equations (E4) and (E5) 


v norm(v=k) 



v norm (v=k) 


sin 0 cos 6 d6 


(E6) 


But 


or 


v = v cos 6 

norm 




v n = V n cos n e 
norm 


(E7) 
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From equations (E6) and (E7), 


/n/2 . 

sin e cos n+ e df) 


n 


= 2v 


n/ -1 „ e n+2A 
( — — cos 6] 

\n + 2 ) 


tt/2 


2 

n + 2 


n 


(E8) 


Now fo r any point on any solid body, with random orientation, the value of 
v norm(v=k) ^ ven by equation (E8) so long as the surface at that point is not shad- 

owed by other parts of the surface of the body. But equation (E8) applies only for a spe- 
cific value of v. 

In the writing of an expression for v ^ orm for a distribution of velocities, use will be 
made of the fact that equation (E8) must apply for each velocity. Use will also be made of 
the statistical frequency of velocity f(v). 

The arithmetic average of the n u power of the normal component of impact velocity, 
then, may be expressed as 


~7t/2 


/ n 

norm 


Jo X 0 d «’ f(v) dv f 0 v £orm(v i k> f < v > dv 




f(v) dv 


r 

Jo 


(E9) 


f(v) dv 


in which the final denominator is the total fractional flux (unity). 
From equations (E8) and (E9) 


or 



(E10) 


v 


n 

norm 


2 

n + 2 


n 


(Ell) 
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APPENDIX F 


AVERAGE VALUE OF A POWER OF IMPACT VELOCITY BY 
INTEGRATION OF BIMODAL LOG-NORMAL EQUATION 


In the main text, expressions were derived for critical values of damage criteria for 
meteoroid impact. Those expressions included an average value of a power of impact 
velocity v°, with n equal to e/3/x. As bimodal log-normal distribution equations have 
been found for velocities Vq and v M , it is therefore desirable to find a solution for the 
function v 11 calculable from those equations. 

The function v n , which involves no negative values of v, may be expressed as 



v n f(v) dv 


For the bimodal log-normal case, equation (3), 



(FI) 


(F2) 


with equal to zero. 

Obviously, equation (F2) can be integrated if the following unimodal expression can 
be integrated 



No simple closed solution or converging series solution has been found for the combina- 
tion of (1) values of n that are not positive integers and (2) nonzero 5. For all other 
cases, however, solutions are found as follows: 

Let 


y = v + 6 


(F4) 


Then, 



or, for positive integral values of n, 



with the conventions that both zero to the zero power and factorial zero equal unity. 
Now, let 


Then 


x = log e y 


(F7) 




(n - i + l)x - — 
2 



dx 


(F8) 


The squared part of the exponent of e may be expanded and terms in x may then be 

2 

combined. The resulting terms in x and x will then form the first two terms of an 
expanded square of a binomial. That square may be completed, so that the exponent of 
e will be a squared binomial and a constant term. The constant power of e may then be 
removed from the integral, with the following result: 



n: 


i! (n - i)! 


(-5) 1 exp i(n 


j(n - i + 1) |m + t 


J (n - i + 1) 


} 
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Now the bracketed part of equation (F10) is an ordinary normal distribution integral 
and may be evaluated as follows: 

Let 


and let 


t = 


x - [p + (n - i + 1)] 
a 


T = 


log e 6 - [p + a 2 (n - i + 1)J 


(Fll) 


(F12) 


Then (see ref. 8 or other textbook on mathematical statistics) 

:k L -( £ ^ l)-t /""'ft*’) 


dt 


= 0. 5 - 


V (-i)i 


2j 


V 2 7 4 -^ 2 j j! (2j + 1) 

1=0 


(F13) 


So, from equations (F10) and (F13), 


__ n , 

= cff V 27 t ^ - — — — - — (-5)* exp|(n - i + 1) |p + -^ a 2 (n - i + 1) 1 

i=0 


x 


QO 

(-i ) 1 — t — 


2 V 2 n / 7 2 j j! (2j + 1) 

j=0 


(6 > 0, n a positive integer) (F14) 


In use of equation (F14) the value of r must be determined for each value of the in- 
dex i with use of equation (F12). The value of the second term within the brackets of 
equation (F14) may then be found in a table of probability integrals (see ref. 8 or other 
textbook on mathematical statistics). 

When 6=0, all parts of the compound series in equation (F14) reduce to zero for 
i > 0. In that case, the value of r from equation (F12) becomes equal to and the 
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bracketed part of equation (F14) becomes equal to unity, so that 


= ccr V 2 H exp|(n + 1) + ^- cr^(n + 1)J j (6 = 0) (F15) 

Equation (F15) is valid for 6=0 and any real value of n. 

If 6 is negative, equation (F12) indicates an imaginary value for r. But, in that 
case, equation (F3) implies nonexistence of values of v less than -6. Hence, the real 
part of the integral in equation (F3) extends from v = -6 to v = °°. Accordingly, the 
real part of the integral in equation (F5) extends from zero to 00 , and the real parts of 
the integrals in equations (F8) to (F10) extend from to +°°. Hence, a value of 
t = -°° may be obtained from equation (F12) and, as in the case with 6 equal to zero 
the bracketed part of equation (F10) becomes equal to unity. If n is a positive integer, 
equation (F10) becomes 


n 


** = coVZrt V — (-6)* expl(n - i + 1) 

Z-J v. {n- i)! I 

i=0 


u 


jju + — cr^(n - i + 1) 


(F16) 


(6 < 0, n a positive integer). 


If 6 * 0, and n is not a positive integer, none of the equations (F14), (F15), or 
(F16) may be used. In that case, equation (F3) must be subjected directly to numerical 
integration. 
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Figure 3. - Average of n* h power of meteoroid impact velocity times concentration factor far from Earth and near Earth. 
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